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Abstract 


It is now easier to simulate catchments and take management measures to reduce flood damage 
with the development and expansion of hydrological tools. the SWAT model was used as a semi- 
distributive model and the SIMHYD model was used as an integrated model to estimate the runoff 
of the Karganrud Talesh catchment. For this purpose, observational data in the period 2006 to 
2018 related to Talesh, Astara, Bandar Anzali synoptic stations and Hashtpar hydrometric station 
were obtained. In modeling by SWAT model, the coefficients N2 and R2 as the objective 
functions in the monthly flow calibration stage were 0.61 and 0.62, respectively, and in the 
validation stage were 0.67 and 0.70, respectively, indicating acceptable results in runoff peak 
simulation and the floods formed in this basin. Acceptable results were also obtained in the 
calibration and validation of the SIMHYD model. The coefficients N2 and R2 in the calibration 
stage were 0.710 and 0.8511, respectively, and in the validation stage were 0.625 and 0.9496, 
respectively. However, due to its semi-distributive nature, the SWAT model has a high accuracy in 
runoff simulation compared to the SIMHYD model, and the outputs obtained from both models 
showed this superiority in flood assessment and estimation, since this model simulates runoff well 
in the daily time period, while the SIMHYD model has satisfactory results in this basin by 
increasing the quality of data in the monthly time period. 
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حکیده 


امروزه با پیشرفت و گسترش ابزارهای هیدرولوژیکی و صرفه‌جویی در زمان» هزینه و 
در شرایط بحرانی به‌منظور اقدامات مدیریتی برای کاهش خسارت‌های ناشی از chee‏ 
استفاده از gee‏ ابزارهایی باعث سهولت در pl‏ شبیه‌سازی حوضه‌های آبریز شده 
است. در پژوهش حاضر از مدل SWAT‏ به‌عنوان یک مدل نیمه‌توزیعی و از مدل 
0 به‌عنوان یک مدل یکیارچه جهت برآورد رواناب حوضه آبریز کرگانرود 
تالش استفاده شده است. بدین منظور از داده‌های مشاهداتی در محدوده زمانی ۲۰۰۹ تا 
۸ مربوط به ایستگاه‌های سینوپتیک تالش, آستارا» بندرانزلی و ایستگاه هیدرومتری 
هشتپر بهره گرفته شده است. در مدل‌سازی توسط مدل SWAT‏ ضرایب N2‏ و R?‏ 
به‌عنوان تابع هدف در مرحله واسنجی دبی ماهانه به ترتیب ۰/۱۱ و ۰/۸۲ و در مرحله 
اعتبارسنجی به ترتیب ۰/0۷ و ۰/۷۰ به دست آمد که نشان‌دهندة نتایج قابل قبولی در 


شببه‌سازی پیک ls)‏ و Bem de‏ شده در این حوضه داشته ینت در واسنجی و 
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اعتبارسنجی مدل SIMHYD‏ نیز نتایج قابل قبولی به دست آمد. به‌گونه‌ای که ضرایب 
2 و R?‏ در مرحله واسنجی به ترتیب ۰/۷۱۰ و ۰/۸۵۱۱ و در مرحله اعتبارسنجی به 
ترتیب ۰/۱۲۵ و ۰/۹۶۹ حاصل شد. مدل SWAT‏ به دلیل ماهیت نیمه‌توزیعی خود از 
دقت بالاتری در شبیه‌سازی سیلاب حوضه موردمطالعه نسبت به مدل SIMHYD‏ 
برخوردار است؛ زیرا این مدل در بازه زمانی روزانه نیز به‌حوبی پیک رواناب را 
شبیه‌سازی می‌کند» درصورتی که مدل SIMHYD‏ با افزايش کیفیت داده‌ها در بازه زمانی 
ماهانه نتایج رضایت بخشی در این حوضه دارد. 

کلیدواژه‌ها: شبیه‌سازی رواناب» مدل SWAT‏ مدل SIMHYD‏ سیل» حوضه 


۱- مقدمه 


برای مدیریت بهتر» اطلاع از بیلان آبی حوضه‌های آبریز ضروری است و به دلیل وقت‌گیر و پرهزینه 
بودن اندازه‌گیری اجزای بیلان آب. مدل‌سازی پدیده‌های هیدرولوژی در حوضه‌آبریز می‌تواند راه‌حل بهینه‌ای 
برای آن‌ها باشد (زارع کاریزی و طالبی. ۱۳۹۵)؛ cpl ple‏ آگاهی دقیق از رفتار هیدرولوژیک حوضه‌های 
آبریز می‌تواند به شبیه‌سازی بهتر این محیط به‌منظور JES‏ مولفه‌های hel‏ بیلان آب کمک شایانی LUS‏ 
(موذن‌زاده و همکاران. ۱۳۹۵). در ASI‏ موارد نحوه توزیع طبیعی جریان‌های سطحی در نواحی مختلف کشور 
با نیازهای زمانی حوضه‌های آبریز منطبق نبوده و ضرورت مهار و کنترل چنین جریان‌های سطحی را به وجود 
آورده است. برآورد حجم رواناب ناشی از نزولات جوی و تغییرات کمی و کیفی آن از مهم‌ترین اقدامات در 
جهت مدیریت و برنامه‌ریزی منابع آب موجود است. امروزه استفاده از روش‌های نوين و غیرمستقیمی 
همچون مدل‌های هیدرولوژیکی برای شبیه‌سازی حوضه‌های آبریز و محاسبه حجم رواناب تولید شده در 
حوضه‌ها به‌عنوان عاملی مهم در وقوع سیل» سبب سهولت در انجام محاسبات و صرفه‌جویی در زمان و 
هزینه شده است. ازاین‌رو کاربرد چنین ابزارهایی در برنامه‌ریزی‌های مدیریتی و اجرای طرح‌های آبی روندی 
اقتانشن ge cad sols‏ ا P‏ کی ان acl i‏ مت اسان یال کور مورت انان 
OLS‏ با توجه به شرایط اقلیمی حاکم از نظر درصد بروز سیل. رتبه بالایی را به خود اختصاص می‌دهد که 
حوضة کرگانرود تالش» واقع در این استان نیز از این yal‏ مستلنی نبوده است؛ زیرا وجود بارندگی‌های منظم 
جوی. مجاورت با دریای خزر و رشته‌کوه‌های تالش باعث تشکیل دامنه‌های جنگلی و پرشیب شده که آن 
منطقه را مستعد پدیدۀ سیلاب کرده است و گسترش روند Of‏ در طول سال‌های متمادی می‌تواند برای 
ساکنین آن منطقه مشکل‌ساز شود. از طرفی دیگر دخل و تصرف غیرمسئولانه و سوء مدیریت در استفاده از 


منابع طبیعی به‌مانند برداشت بی‌روية درختان جنگلی» تغییرات کاربری اراضی و احداث سازه‌های انسانی در 
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بستر و حریم رودخانه. در بروز رواناب مازاد بر نفوذ و درنتیجه مخاطرةٌ سیل بسیار مژثر بوده است. 
شهرستان تالش از نظر موقعیت جغرافیایی خود در بین دو واحد بزرگ ساختمانی البرز-تالش و دریای خزره 
bow‏ مناسبی را برای تشکیل اشکال و پدیده‌های ژئومورفولوژیکی فراهم ساخته است. بخصوص با دارا 
بودن شرایط اقلیمی و محیطی متفاوت در طول سال‌های گذشته دچار خسارات فراوان از آسیب‌های طبیعی 
گردیده که بدون شک پدیدة سیل به‌عنوان یکی از جدی‌ترین خطرات و بلایای طبیعی این منطقه محسوب 
می‌شود. در طول سالیان متمادی مردم این نواحی به دلیل ویژگی‌هایی همچون: مرغوبیت خاک منطقه و 
ذخیره آب کافی» اغلب بدون در نظر گرفتن دورة بازگشت و حجم رواناب تولید شده در هرسال و با وجود 
تهدیدهای حاکم» در hilo‏ رودخانه‌ها سکونت گزیده و گسترش یافته‌اند؛ اما با توجه به اهمیت 
رودخانه‌های استان OLS‏ و doe‏ اقتصادی فراوان daol‏ تاکنون مطالعات ویژه‌ای بر روی آن‌ها انجام نشده و 
متأسفانه در هيچ‌یک از پژوهش‌های انجام شده ساماندهی و مدیریت سیلاب به‌صورت جدی صورت نگرفته 
cul‏ مدل SWAT‏ یکی از ابزارهای تحلیلی برای شبیه‌سازی فرآیندهای هیدرولوژیکی در حوضه‌های آبریز 
پیچیده است که با توجه به تغییرات خاک کاربری اراضی» کیفیت آب و شرایط اقلیمی در طولانی‌مدت و 
حوضه‌هایی که فاقد آماربرداری منظم هستند. کاربرد دارد (گاسمن و همکاران؛ ۲۰۰۷). مدل SIMHYD‏ نیز 
به‌عنوان مدل مفهومی ساده بارش -رواناب به شکلی یکپارچه عمل کرده و در مطالعات مربوط به شبیه‌سازی 
رواناب (پیل و همکاران ۲۰۰۲ usb‏ تغییرات اقلیمی بر رواناب (شن و همکاران؛ ۲۰۰۷) و ارزیابی 
منطقه‌ای ger)‏ و سیریواردنا ۲۰۰۵) به‌صورت موفقیت‌آمیزی عمل کرده است. امروزه این ابزارها بسیار 
موردتوجه جغرافیدانان و مهندسین منابع آب Jel‏ و خارجی در بحث مدیریت منابع آب واقع شده است. 
حبیبی و گودرزی (۱۳۹۷) در تحقیق خود کاربرد مدل نیمه توزیعی SWAT‏ در شبیه‌سازی رواناب 
حوضه حبله‌رود را بررسی کرده‌اند. بررسی‌های آن‌ها OUS‏ داد که از مدل SWAT‏ می‌توان برای شبیه‌سازی 
رواناب در مقیاس زمانی ماهانه در حوضة موردنظر برای برنامه‌ریزی‌های منابع آبی و اقدامات مدیریتی لازم 
در این حوضه به‌خوبی استفاده نمود. دستحردی و همکاران (۱۳۹۸) در پژوهش خود کارایی دو مدل 
هیدرولوژیکی AWBM‏ و 511۳17۳ و شبکه عصبی MID)‏ و (Rbf‏ را برای شبیه‌سازی بارش-رواناب در 
منطقه حوضه باراریه نیشابور مقایسه کرده و GLE‏ حاصل OLS‏ می‌دهد که تطابق زیادی بین مقادیر اوج 
بارندگی با اوج رواناب شبیه‌سازی شده توسط مدل‌ها و رواناب حداکثر خحروجی از حوضه. وجود دارد. 
گودرزی و فاتحی‌فر (YIA)‏ با استفاده از مدل هیدرولوژیکی SWAT‏ در GIS lcs‏ در پهنه‌بندی خطر 
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سیلاب در اثر تغییرات اقلیمی تحت سناریو 8.5 RCP‏ مطالعه‌ای در Lege‏ آذرشهر چای انجام داده‌اند و 
نتایج حاصل از تحقیق آن‌ها حاکی از آن است که ممکن است در یک زیرحوضة کوچک در اواسط حوضه به 
جهت داشتن شرایط مستعد سیل Jub‏ شیب فراوان داشتن شکلی گرد. خاک نفوذناپذیر یا پوشش‌گیاهی 
ضعیف. در اثر باران‌های شدید اقلیمی دبی سیلابی Soy‏ به وقوع بپیوندد. قضاوی و همکاران(۱۳۹۹) در 
ub‏ پژوهشی به بررسی کارایی مدل SWAT‏ در برآورد دبی‌های روزانه حوضه‌های فاقد آمار با رویکرد 
منطقه‌بندی در مناطق خشک پرداخته و نتایج حاصل از پژوهش OUS OLT‏ می‌دهد که خروجی‌های به دست 
آمده موفقیت‌آمیز و قابل‌قبول بوده و استخراج این نتایج می‌تواند به درک و شناحت بهتری از رفتار 
هیدرولوژیکی حوضه‌های BE‏ آماربرداری برای اهدافی چون برنامه‌ریزی و مدیریت منابع آب ازجمله کنترل 
سیل و تأمین آب شرب و کشاورزی منجر شود. بالویی و همکاران(۱2۰۰) کارایی مدل SWAT‏ را در 
شبیه‌سازی رواناب و رطوبت موجود در خاک در حوضه دویرح استان pL!‏ بررسی نموده‌اند. این مدل در 
شبیه‌سازی Oy‏ آبی موجود در این حوضه عملکرد بسیار موثر داشته و قابل استفاده برای پردازش phe‏ 
ویژگی‌های حوضه دویرج نیز می‌باشد. اسکوبار و همکاران" (۲۰۱۵) حوضه آبریز مولاوین» به‌عنوان 
ذخیره‌گاه جنگلی KASU‏ را با استفاده از پنج مدل مفهومی توده‌ای بارش-رواناب شبیه‌سازی کرده و آن‌ها 
پس از مدل‌سازی با توجه به نتایج خود پنج الگوریتم موجود در نرم‌افزار RRL‏ را از بهترین تا ضعیف‌ترین 
طبقه‌بندی نمودند که آرایش مدل‌ها عبارتند از: AWBM ; SMAR SACRAMENTO SIMHYD Tank‏ 
هونگیان و همکاران" (۲۰۱۷) به ارزیابی در رابطه با توانایی و سازگاری مدل SIMHYD‏ در حوضه ;4-55 
سونگوا پرداعتند. آن‌ها دریافتند که استفاده از این مدل در مناطق نیمه‌مرطوب te‏ از مناطق نیمه‌عشک است. 
این مدل همچنین دو مکانیسم تولید رواناب » ذخیره‌سازی و نفوذ بیش‌ازحد را در نظر می‌گیرد . با ساختار 
ساده و محاسبه راحت نیاز به پارامترهای کمتری دارد. اما این مدل در اوج سیل دارای یک خطای شبیه‌سازی 
بزرگ است. با توجه به تأثیر توزیع نابرابر باران و جاری شدن سیل» حوضه به چندین بلوک تقسیم می‌شود. 
درنتیجه هر بلوک برای محاسبه تلاقی به‌طور جداگانه به چندین واحد تقسیم شده و سرانجام به‌صورت خطی 
در ایستگاه خروجی قرار می گیرند که ازاین‌جهت می‌تواند صحت مدل را بیشتر ope‏ ببخشد و از این نظر 
برای ارتقا بیشتر مدل در حوضه رودخانه سونگوا از اهمیت بسیاری برخوردار است. یو و همکاران" (۲۰۱۸) 
مقاله خود را با عنوان بهبود مدل SWAT‏ برای شبیه‌سازی مبتنی بر رویداد سیل در بازه زمانی روزانه ارائه 
نمودند. نتایج درجه‌بندی و اعتبارسنجی داده‌های روزانه نشان‌دهندة بهبود عملکرد شبیه‌سازی به‌ویژه 
شبیه‌سازی پیک سیلاب است. آن‌ها همچنین معتقدند بهبود این مدل» پایه و اساس رسیدگی به مسائل مربوط 
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به کیفیت CT‏ ناشی از وقوع سیل خواهد بود. ol‏ و عابد" (۲۰۲۱) رواناب سطحی در بحرالنجف را با 
استفاده از مدل SWAT‏ برآورد نموده‌اند. سه سناریو فرضی بر اساس میزان بارندگی مفروض اجرا شد که 
نتایج نشان دادند سه سطح Ol‏ منجر به پرشدن دریاچه محصور کل منطقه موردمطالعه و فراتر رفتن از مرزها 
و نیز باعث سیل بخشی از شهر باستانی نجف شده است. سیلوا دیاس و همکاران" (۲۰۲۱) عملکرد مدل‌های 
هیدرولوژیکی در شبیه‌سازی جریان رودخانه‌ای را با استفاده از مدل‌های SWAT‏ و HEC-HMS‏ ارزیابی 
کرده و به این نتیجه رسیدند که با توجه به شاحص‌های به دست edel‏ هر دو سیستم نرم‌افزاری به‌طور 
رضایت‌بخشی قادر به شبیه‌سازی جریان بودند. باسو و همکاران" (۲۰۲۱)» ابوزریگ و بنی هانی * (۲۰۲۱) و 
دالیسون و همکاران" (۲۰۲۲) نیز با استفاده از مدل SWAT‏ در پژوهش خود توانستند به نتایج خوبی دست 
dub‏ تحت شرایط این SWAT «JU.‏ نسبت به HEC-HMS‏ برتری دارد و از JE‏ آماری برای 


پیش‌بینی‌های آینده مناسب‌تر است. 
-Y‏ مواد و روش 
۱-۲- موقعیت جغرافیایی حوضه آبریز کرگانرود 


حوضه آبریز کرگانرود با مساحتی بالغ بر 1۱٥/٤‏ کیلومترمربع و با محیطی حدود ۱۲۷ کیلومتر یکی از 
رودخانه‌های شهرستان تالش و همین‌طور از مهم‌ترین و پرآب‌ترین رودخانه‌های استان گیلان می‌باشد که در 
شمال‌غرب این استان از ارتفاعات ۲۰۰۰ تا ۲۵۰۰ متری تالش با مختصات جغرافیایی YVD ۶۱ ٣٤‏ تا VE‏ 


YVL] ۸‏ عرض شمالی ۶۸5 ۳۱ EAL‏ تا EAL ۵۸ £A‏ طول‌شرقی سرجشمه گرفته است. از نظر تقسیمات 
سیاسی این حوضه از شمال به حوضه هره‌دشت لیسار» از جنوب به حوضه ناورود اسالم از غرب به استان 


اردبیل حوضه آبریز آریاچای و هروآباد و از شرق به دریای خزر محدود می‌شود (شکل AV‏ 


1 Farhan and Abed 

2 Silva Dias et al 

3 Basu et al 

4 Abu-Zreig and Bani Hani 
5 Dallison et al 


YA‏ جغرافیا و مخاطرات محیطی شمارة اول 


شکل ۱- نقشه موقعیت جغرافیایی حوضه آبریز گانرود تالش 


۲-۲- معرفی مدل SWAT‏ 


۲ مخفف عبارت Soil and Water Assessment Tools‏ می‌باشد که یک مدل پیوسته و نیمه 
توزیعی بوده و توسط جف آرنولد (آرنولد و همکاران» OA‏ برای سرویس تحقیقاتی کشاورزی آمریکا 
توسعه یافته است. این مدل برای ارزیابی اثرات کوتاه‌مدت و بلندمدت مدیریت بر روی آب» رسوب. مواد 
شیمیایی و کشاورزی در حوضه‌ها انجام می‌شود. مدل SWAT‏ دارای cle‏ فیزیکی بوده و شامل دو بخش 
اصلی است. در بخش اول میزان آب» رسوب. فرسایش خاک رشد گیاه» مواد مغذی و آفت‌کش‌هایی که وارد 
al yl‏ اصلی می‌شوند را کنترل می‌کند و در بخش دوم روندیابی و جهت جریان آب. حرکت رسوبات و 
شرایط آب و هوایی در یک حوضۀ آبریز شبیه‌سازی می‌شود. تلفیق GIS‏ با مدل نیمه فیزیکی SWAT‏ باعث 
شده است تا قدرت تحلیل لایه‌های اطلاعاتی به‌منظور ایجاد لایه‌های جدید افزایش LE‏ و در مدیریت 
سیستم‌های پیچیده آبی نقشی کاربردی ایفا کند. چرخه هیدرولوژیکی این مدل بر معادله بیلان آبی هر HRU‏ 
استوار است (آبابایی و سهرابی. ۱۳۸۸). به‌طورکلی معادله بیلان آبی» اصلی‌ترین معادله در هر نوع فرآیند 


SW, = SW, + CT m Qsurf m Eq m Weep a Qow) 


1 Arnold et al 
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در این SW, Soles‏ ذخیره رطوبت نهایی در خاک t (mm)‏ زمان SWo‏ مقدار adsl‏ آب در خاک mm)‏ 
Raay‏ ميزان بارش روزانه» Ol x Qsurf‏ روانات سطحی روزانه ) Ol x Ea «(mm‏ تبخیر 3 تعرق روزانه (mm)‏ 
Ope Was‏ آبی که از پروفیل خاک به منطقه غیراشباع نفوذ می‌کند (mm)‏ و ,یم میزان جریان نفوذی در 


سفره زیرزمینی به‌صورت روزانه «(mm)‏ می‌باشد (نیچ و همکاران» 0+ (Ye‏ 
-Y-Y‏ داده‌های ورودی مدل 


ورودی‌های مدل SWAT‏ برای انجام شبیه‌سازی شامل نقشه رقومی ارتفاع (DEM)‏ با قدرت تفکیک 
کے ۱۲/۵ at ee‏ شیب cle UU Sb‏ او a BCS hoy‏ شاک برش cida‏ 
shes‏ کاک Ohi‏ وه was‏ اش و او alat‏ شا سای ارام ada‏ تست 
مربوط به dle‏ ۲۰۱۸ می‌باشد که با استفاده از روش تصحیح اتمسفری سریع (QUAC)‏ و روش حذاکثر 
احتمال در la‏ ۷1 ترسیم شده و با نتایج ضریب (Kappa= 0/97) LE‏ و دقت Overal ) AS‏ 
(accuracy= 0/98‏ بالا به دست آمده است (شکل ۲). همچنین از آمار و اطلاعات روزانه و ilal‏ ماکزیمم و 
مینیمم دما بارش» رطوبت نسبی. سرعت با تابش خورشیدی مربوط به سه ایستگاه تالش آستارا و 


بندرانزلی به‌عنوان ورودی‌های این مدل نیز استفاده شده است (حدول 6 
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شکل (A -Y‏ نقشه شیب. 8) کاربری اراضی و (C‏ خاک حوضه کرگانرود به‌عنوان ورودی مدل SWAT‏ 


1 Neitch et al 


۳۰ جغرافیا و مخاطرات محیطی شمارة اول 


حدول ۱- اطلاعات ایستگاه‌های منتخب برای مدل‌سازی 


نام ایستگاه طول جغرافیایی | عرض جغرافیایی | ارتفاع از سطح دریا | نوع ایستگاه | دوره آماری 


-Yen سینوپتیک‎ ۷ ۳۷۲۱۲۱۵۰۲۸ | £ALI£a'oY"E تالش‎ 
-ye بندر آستارا | ۸۱۵۱۱۵۲۴8 | ۳۸۲۱۵۶۳۲۱۳۲ ۷ سینوپتیک‎ 
-Ye سینوپتیک‎ -WwA YVD £VNA"N | 1٩1۲۷۲۷۳۳ | Jl بندر‎ 
-Yed هیدرومتری‎ ۳۵ YVO ۶۹۸۲۸ 0۵ هس‎ 


در این GES‏ برای بررسی میزان سیل‌خیزی حوضه آبریز کرگانرود از پارامترهای شماره منحنی (CN)‏ در 
بازه صفر تا ۰۰ Qiu‏ ميزان رواناب سطحی روزانه Raay (mm)‏ ميزان بارش روزانه و 5 پارامتر نگهداشت 
رطوبت (mm)‏ مطابق بر معادله زیر استفاده شده است: 


s = 25/5 ŒŒ — 10) 


SIMHYD مدل‎ -£-Y 


GRRE" jailed‏ استفاده از داده‌های مشاهداتی بارش تبخیر و تعر و دبی مشاهناتی برای شیه‌سازی 
یک مدل بارش -رواناب تهیه و تدوین شده است. این نرم‌افزار حاوی پنج مدل بارش -رواناب توده‌ای «gale‏ 
A A VEDE EE SN En EEE Ng A E‏ 
SMAR TANK SIMHYD 7‏ و SACRAMENTA‏ مدل مفهومی یکپار i23 SIMHYD' im‏ 
ساده‌شدة نرم‌افزار HYDROLOG‏ مربوط به سال ۱۹۷۲ و نرم‌افزار MODHYDROLOG‏ مربوط به سال 
1a‏ اش یکی adits See a‏ با تسس ماش سر نیتم NOB‏ ای سل تاه 
ارزیابی رواناب از سه منبع استفاده می کند که عبارتند از: رواناب مازاد بر نفوذه جریان زیرسطحی و رواناب 
مازاد بر اشباع و جریان پایه. میزان تغذیه آب زیرزمینی که تابعی خطی از رطوبت خاک می‌باشد نیز توسط 
مدل SIMHYD‏ برآورد می‌گردد (پودگر" ۲۰۰۳). درواقم عملکرد مدل SIMHYD‏ به این صورت است که 
ابتدا ذخایر خالی‌شدۀ برگاب» توسط تبخیروتعرق به‌وسیلهٌ بارش روزانه پر می‌شود. سپس بارش مازاد تحت 
puc‏ یک تابع نفوذ قرارگرفته و ظرفیت نفوذ را مشخص می‌نماید. بارش مازادی که از ظرفیت نفوذ فراتر رفته 
و تجاوز می‌کند. رواناب را به وجود می‌آورد. رطوبت نفوذ کرده تحت تأثیر یک تابع رطوبتی خاک قرارگرفته 
که این رطوبت به جریان‌های داخلی. آب‌های زیرزمینی و ذخایر رطوبت خاک می‌پیوندد؛ cpl ple‏ این مدل 


بر اساس داده‌های مشاهداتی بارش تبخیروتعرق 3 دبی رودخانه» می تواند رواناب را در محد‌وده زمانی 


1 Rainfall Runoff Library 
2 Simple Hydrology 
3 Geoff Podger 
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موردنظر شبیه‌سازی کند که در این تحقیق از داده‌های مشاهداتی به‌صورت روزانه در $56 زمانی ۱۳ ساله بهره 


۵-۲- داده‌های ورودی مدل SIMHYD‏ 


ازآنجایی که مدل SIMHYD‏ نسبت به مدل SWAT‏ از ورودی‌های کمتری برخوردار است براین اساس 
سه نوع ورودی ازجمله داده‌های تبخیروتعرق و میزان بارش از سازمان CT‏ و هواشناسی استان گیلان و 
داده‌های مشاهداتی مربوط به دبی جریان از سازمان آب‌منطقه‌ای این استان مربوط به ایستگاه سینوپتیک تالش 
و ایستگاه هیدرومتری هشتیر که در قسمت خروجی رودخانه قرارگرفته است. در بازه زمانی ۲۰۰۲ تا ۲۰۱۸ 
تهیه گردید و به فرمت مخصوص نرم‌افزار مرتب شده. سپس برای مدل تعریف شده است. به دلیل تازه 
تأسیس بودن ایستگاه مورد استفاده در شبیه‌سازی و همچنین عدم وجود سنوات آماری کافی بلندمدت. لذا 
این محدوده آماری برای انجام مدل‌سازی در نظر گرفته شده است. در این مدل ضریب نش- ساتکلیف 
(EN2)‏ به‌عنوان تابع هدف که اختلاف Ole‏ داده‌های مشاهداتی با مقادیر شبیه‌سازی شده را نمایش می‌دهد 
ضریب تبیین RB)‏ و ريشه میانگین مربعات خطا (RMSE)‏ نیز برای مراحل واسنجی و اعتبارسنجی در نظر 


گرفته شد. 
-Y‏ نتایج و Cou‏ 
۱-۳- اجرای مدل و شبیه‌سازی هیدرولوژیکی رواناب توسط SWAT-CUP‏ 


پس از تعریف کردن تمامی نقشه‌ها و داده‌های ورودی به مدل. همچنین در صورت SLE‏ برای دریافت 
نتایجی با دقت بالا مدیریت‌های مربوط به SE‏ کاربری اراضی, آبیاری يا در صورت وجود سد و غیره 
اطلاعات مربوط به آن در مدل اعمال شده سپس در مرحله نهایی تاریخ شروع و اتمام دوره شبیه‌سازی از 
سال ۲۰۰۱ تا ۲۰۱۸ تعداد سال‌های گرم شدن مدل" به مدت سه سال تعیین گردید و مرحله مدل‌سازی 
به‌صورت روزانه اجرا شد. یکسری اطلاعات خروجی شبیه‌سازی شده به شکل فایل Textinout‏ تولیده ott‏ 
که این خروجی‌ها به نرم‌افزار SWAT-CUP‏ به‌منظور انجام فرآیند واسنجی و اعتبارسنجی با تعیین 
الگوریتمی که مدل بر اساس آن کالیبره خواهد شد افزوده شدند. با توجه به مدل‌سازی انجام گرفته و نتایج 
حاصل از گزارش مدل SWAT‏ میانگین تبخیروتعرق را به میزان ٤1٩‏ میلی‌متر: میانگین بارش ۹٩۱۲‏ میلی‌متر 
در حوضه آبریز کرگانرود محاسبه نموده است. با در نظر گرفتن اقلیم حاکم بر حوضه کرگانرود تالش. نوع 


1 Warm Up 


۳ جغرافیا و مخاطرات محیطی Silas‏ اول 


کشت برنج در بخش کشاورزی کرگانرود به‌عنوان گیاهی با نیاز آبی بسیار بالا و بیشترین تغذیه توسط 
بارندگی و درنهایت با توجه به گروه هیدرولوژیکی خاک‌های واقم در حوضه. این‌طور استنباط می‌شود که 


polis‏ به دست آمده با مقادیر موجود در واقعیت هماهنگی داشته و منطقی به نظر می‌رسد. 
۲-۳- واسنجی و اعتبارسنحی مدل SWAT‏ 


جهت ope‏ نتایج شبیه‌سازی و اطمینان دقیق از شباهت و هماهنگی بالای مقادیر بارش -رواناب تولید 
شده با مقادیر موجود در واقعیت ضرورت انجام واسنجی و اعتبارسنجی مدل شبیه‌سازی شده فراهم می‌شود. 
برای این کار در ابتدا پس از تعریف 2 655 Textinout‏ به‌عنوان خحروجی حاصل از مدل‌سازی در SWAT‏ در 
محیط SWAT-CUP‏ افزوده شده و الگوریتم SUFI‏ به‌عنوان روش منتخب برای انجام کالیبراسیون انتخاب 
شده است؛ زیرا این روش به‌راحتی می‌تواند حجم وسیعی از داده‌های مشاهداتی را پردازش کرده و امکان 
کالیبراسپون همزمان پارامترها بر ah‏ داده‌های ارزیابی شده در متغیرهای متفاوت را دارا می‌باشد. سپس بعد از 
وارد کردن داده‌های دبی رودخانه با فرمت نرم‌افزان در ابتدا پارامترهای احتمالی با توجه به بررسی تجربیات 
پژوهش‌های گذشته توسط دیگران مانند بویتیاس و همکاران (۲۰۱۷)» شکوهی‌فر و همکاران: (۲۰۲۱): 
نیکخو امیری و همکاران(۱۳۹۸). طالبی و همکاران(۱۳۹۹) جهت برآورد = OL‏ جمع‌آوری شده سپس به 
سوات کاپ معرفی و آنالیز حساسیت برای ۲۵۵۹ داده ب‌صورت روزانه انجام گرفت. دامنه‌ای از polis‏ 
حداکثر و حداقل برای هر پارامتر در نظر گرفته شد و سپس با بررسی خروجی کالیبراسیون و تولید ضریب 
نش مناسب. از میان تمام پارامترهای آنالیز شده ۲۳ پارامتر موثر با حساسیت بالا در تولید رواناب با ۵۰۰ بار 
تکرار در شبیه‌سازی و نمونه‌برداری تعیین گردید و تثبیت شد. این پارامترهای موثر در تولید رواناب در 


(Y J te)‏ نمایش داده شده است. 


1 Boithias et al 
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جدول -Y‏ پارامترهای مؤثر در تولید رواناب 


درجه حساسیت نام پارامتر شرح پارامتر s‏ حداقل 
CN2 \‏ شماره منحنی ۲ | -A‏ 
CH_N2 Y‏ ضریب مانینگ در رودخانه اصلی TAI,‏ ۰ 
CH-K2 Y‏ ضریب_هدایت هیدرولیکی مؤثر کف رودخانه اصلی ۱۳۰ As‏ 
SFTMP 1‏ متوسط دمای هوا برای bas‏ _باران به برف (C^)‏ ۰ 0- 
TIMP ۵‏ فاکتور تأخیر دمای توده برف ۵ | ۰/۲۵ 
SOL_K 1‏ ضریب _هدایت هیدرولیکی خاک =/A =X (mm/hr)‏ 
GWQMN ۷‏ لایه یا برگشتی اتفاق بیفتد (MHO‏ کم عمق برای این که ۰ | ۲۸۰۰ 
SOL_AWC A‏ ظرفیت cl‏ قابل دسترس خاک —+/Y (mmH2O/mm)‏ ۰/۷ 
SOL_BD q‏ جرم مخصوص خاک( =A | ۲ (gtlom?‏ 
RCHRG_DP E‏ ضریب نفوذ آبخوان ۵ | ۰/۲۵ 
CODICEM ETT ALPHA BNK M‏ \ 0/ 
ALPHA_BF 1۲‏ جریان یایه 8 ضریب \ 0/+ 
SMFMN ۳‏ نرخ ذوب برف در ۲۱ دسامبر E (mmH20/C?- day)‏ ۵ 
GW_REVAP H‏ ضریب برگشت آب زیرزمینی از AY‏ کم عمق به لایه UN‏ ۰/۱ 
CANMX ۱۵‏ حداکثر آب ذخیره شده در تاج گیاه برای مرتم Os‏ ۰ 
SLSUBBSN M‏ متوسط طول شیب =A Eg (m)‏ 
ESCO ۷‏ ضریب جبران تبخیرخاک \ UN‏ 
CANMX \A‏ حداکثر آب ذخیره شده در تاج گیاه برای برنج *0 . 
CANMX E‏ حداکثر آب ذخیره شده در تاج گیاه برای جنگل *0 . 
DEEPST Ye‏ عمق آب در آبخوان عمیق موی | Vers‏ 
SMFMX YA‏ نرخ ذوب برف در ۲۱ ژوئن (mmH20O/C*- day)‏ ۵ . 
SMTMP YY‏ دمای ذوب توده برف (C°)‏ ۵ ۰ 
EPCO yy‏ ضریب جبران جذب گیاه \ UN‏ 


به‌طورکلی یک دورۀ ۱۳ dle‏ برای شبیه‌سازی رواناب به شکل روزانه تعیین گردید. بدین‌گونه که 
داده‌های مربوط به سال ۲۰۰۳ تا ۲۰۰۸ برای دوره سازگار شدن مدل» از سال ۲۰۰۹ تا ۲۰۱۵ یک دوره آماری 
۷ساله برای اجرای فرآیند واسنجی مدل و یک دوره سه‌ساله از سال Y YA‏ تا ۲۰۱۸ برای فرآیند اعتبارسنجی 
در نظر گرفته شد. همچنین در این فرآیند ضریب نش ساتکلیف(5 به‌عنوان تابع هدف در نظر گرفته شده 
است که میزان اختلاف bly,‏ مشاهداتی با رواناب شبیه‌سازی شده را OUS‏ می‌دهد. 


۳-۳- نتایج حاصل از انجام فرآیند شبیه‌سازی رواناب روزانه و ماهانه مدل SWAT‏ 


در فرآیند واسنجی تغییر مقادیر پارامترها Lb‏ به‌گونه‌ای انجام شود که در نمودار شبیه‌سازی شده علاوه 
بر تولید —- نش ساتکلیف مناسب. نتایج تا حد زیادی در محدوده باند عدم قطعیت .407 واقع شود. به 
همین دلیل هر چه ضریب p factor‏ 4 یک و ضریب ar factor‏ صفر نزدیک‌تر باشد نتایج Baal cass ud‏ 


از دقت بالایی در شبیه‌سازی برخوردار است. اجرای واسنجی تا زمانی ادامه می‌یابد که بر اساس الگوریتم 


Yt‏ جغرافیا و مخاطرات محیطی شمارة اول 


انتخاب شده ضرایب موردنیاز با مقادیر قابل قبولی محاسبه شوند. در (جدول و شکل (Y‏ نتایج حاصل 5l‏ 


فرآیند کالیبراسیون نمایش oala‏ شده است. 


جدول ۳- نتایج دقت ضرایب ارزیابی مدل در مرحله واسنجی دبی روزانه 


نتایج شاخص‌های ارزیابی مدل | ضریب نش- ساتکلیف(8) | ضریب Roped‏ ضر یب r-factor‏ ضریب p_factor‏ 
ZY ۱/۳۰ ۲ ۰/0۰‏ 


مرحله واسنجی 


mum ۲ 
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شکل ۳- نمودار باند عدم قطعیت برای مرحله واسنحی دبی روزانه o>‏ ضه کرگانرود 


پس‌ازاینکه بر روی مدل شبیه‌سازی شده فرآیند کالیبراسیون انجام گرفته و پارامترها تعدیل گردید. مدل 
در بازه زمانی مستقلی از مر al>‏ واسنجی. L‏ داده‌های مشاهداتی مربوط به سال‌های LI>‏ به فرمت نرم‌افزار 
تعریف شده و فرآیند اعتبارسنجی انجام II‏ در این al> „o‏ داده‌های دبی روزانه حوضه کرگانرود از سال 
٩‏ تا ۲۰۱۸ مربوط به ایستگاه هیدرومتری هشتیر با همان ۲۳ پارامتر مرحله واسنجی و در همان بازۀ 
مقادیر حداکثر و حداقل با تکرار ۰ بار در هر شبیه‌سازی استفاده شده است. نتایج این al>‏ در (حدول و 
شکل ۶) نشان داده شده است. 


حدول e5-t‏ ضرایب به دست آمده در مرحله اعتبار سنجی دبی روزانه 


نتایج شاخص‌های ارزیابی مدل ضریب نش - (NS)CAJSSL.‏ ضریب Roes‏ | ضریب ۲-۵٤0۲‏ ضریب p_factor‏ 
مرحله اعتبارسنجی ۰/۳ ۰/۹۰ AR!‏ ۰/۷۱ 
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شکل ۶- نمودار باند عدم قطعیت برای مرحله اعتبارسنجی دبی روزانه حوضه کرگانرود 


همان‌گونه که از نمودارهای خروجی حاصل از واسنجی و اعتبارسنجی استنباط می‌شود. مدل polis‏ 
شبیه‌سازی را در برخی روزها بیشتر از polis‏ مشاهداتی برآورد کرده است و polio‏ شبیه‌سازی شده از حالت 
یکنواحتی پیروی می‌کنند. از علت‌های اصلی این موضوع از نگاه شن و همکاران" »)۲۰۱٤(‏ وجود محدودیت 
در مدل SWAT‏ در توزیع بارش از لحاظ مکانی است. به‌طور مثال اگر در بخشی از حوضۂ آبریز بارندگی با 
شدت بیشتری رخ دهد مدل این حجم از بارش را در تمام زیرحوضه‌ها JU‏ می‌دهد و درنهایت باعث 
شبیه‌سازی دبی بیش از مقادیر مشاهداتی می‌شود. از دیگر عوامل تأثیرگذار براین موضوع عدم وجود 
ایستگاه‌های کافی در منطقه است که عواملی چون توپوگرافی‌های متنوع ازجمله کوهستانی و جنگلی بودن 
منطقه. عدم نیاز سازمان‌های هواشناسی و آب منطقه‌ای به ایستگاه‌های بیشتر از دلایل کمبود ایستگاه‌های 
موردنیاز برای مدل است. از منظری دیگر کیفیت پایین داده‌های ثبت شده توسط سازمان‌ها يا وجود داده‌های 
مفقود شدة بسیاری در میان ورودی‌ها از مهم‌ترین دلایل مژثر در خروجی مدل است که مانع از شبیه‌سازی 
iso‏ مطابق با واقعیت می‌شود. برای اطمینان از نتایج به دست آمده و همچنین dallas‏ پژوهش‌های مشابه 
انجام شده در فرآیند واسنجی و اعتبارسنجی به‌صورت روزانه» بار So‏ داده‌های مشاهداتی به شکل ماهانه 
تعریف شده و Glee‏ با شبیه‌سازی روزانه از همان پارامترها با مقادیر تعیین شده استفاده شده است؛ زیرا 
وجود تعداد بسیار زیاد و متنوع داده‌های مشاهداتی در مدل‌سازی روزانه. کیفیت Gal‏ و يا وجود مقادیر 
مفقود شده در داده‌های روزانه ثبت شده توسط ایستگاه‌ها از دلایلی بشمار July‏ که می‌توانند نتایج 
مدل‌سازی را تحت تأثیر خود قرار دهند؛ بنابراین فرآیند واسنجی بار دیگر در محدوده زمانی ۲۰۰۳ تا ۲۰۱۵ 
با در نظر گرفتن سه سال اول برای سازگار شدن مدل (Warmup)‏ و سه سال آخر نیز از سال ۲۰۰۱ تا 


۸ برای اعتبار نیح به‌صورت ماهانه انجام شده است (جدول و شکل ۵ و AV‏ 


1 Chen et al 


۳۹ جغرافیا و مخاطرات محیطی شمارة اول 


حدول 0— نتایج دقت ضرایب ارزیابی ماهائه مدل در مرحله واسنحی oe‏ ماهانه 


نتایج شاخص‌های ارزیاپی مدل | ضریب نش- (NS) AAS be‏ ضریب Que‏ ضریب r-factor‏ ضریب p. factor‏ 
مرحله واسنجی Y UM‏ \0/\ ۰/۹ 
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شکل ۵- نمودار AL‏ عدم قطعیت ماهانه برای مرحله واسنحی eo‏ ماهانه حور ضه کر گانرود 


جدول 7 نتایج ضرایب ماهانه به دست آمده در مرحله اعتبارسنجی دبی ماهانه 


p. factor ضریب‎ r-factor ضر یب‎ Queue نتایج شاخص‌های ارزیابی مدل | ضریب نش- ساتکلیف(5) ضریب‎ 
YAN \/£Y ۰/۸۷۰ JAVA مرحله اعتبارسنجی‎ 
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شکل -A‏ نمودار باند عدم قطعیت ماهانه برای مرحله اعتبارسنجی دبی ماهانه حوضه کر گانرود 


با نگاه اجمالی به مقادیر به دست آمده در شبیه‌سازی رواناب به‌صورت روزانه و ماهانه مطابق با نظریه 
موریاسی و همکاران" Qe)‏ عباسیور و همکاران QV)‏ آی ال" (۲۰۱۷) و کرمی و بیاتی 


خطیبی(۱۳۹۸). E‏ مفدار تابع هدف در نظر گرفته شده به‌عبارت‌دیگر ضریب نش ساتکلیف(N2)‏ ميان 


1 Moriasi et al 
2 Ayele et al 
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--۱/۵ ارزیابی شود نتایج به دست skal‏ رضایت بخش خواهد بود و اگر مقدار تابع هدف ميان (Ago SI‏ 
باشد در ردیف نتایج خوب قرار می‌گیرد؛ بنابراین با توجه به این موضوع poli‏ حاصل از فرآیند واسنجی و 


اعتبارسنجی در این مطالعه موردپذیرش بوده و از نتایج خوبی برخوردار است (جدول AV‏ 


جدول ۷- مقادیر و طبقات ضرایب R‏ و NS‏ 


طبقات polis‏ | شاخص طبقات PE‏ شاخص 
ناکارآمد «Ji‏ ناکارآمد ۰/۵ 
قابل‌قبول 0 — رضایت‌بخش | ۰/٩‏ — ۰/۵ 
رضایت‌بخش | NS vio‏ خوب ۷ ۰۱/۱۰ R?‏ 
US‏ ۰/۸۷۵ 
= خیلی خوب | ۰/۷۱ 
خیلی خوب | ۱ - 


منبع: موریاسی و همکاران(۲۰۰۷)» عباسپور و همکاران(۲۰۰۷)» آی آل(۲۰۱۷) و کرمی و بیاتی خطیبی(۱۳۹۸) 
1-۳ - واسنجی و اعتبارسنحی مدل SIMHYD‏ 


تحلیل حساسیت مدل نسبت به یک پارامتر خاص ازجمله عوامل بسیار مهم و تأثیرگذار در روند 
مدل‌سازی است. چراکه از این طریق می‌توان میزان تأثیر یک پارامتر را نسبت به دیگر پارامترها مشخص نمود 
و بهترین خروجی را به دست آورد؛ بنابراین پارامترهایی با حساسیت بالا بعد از انجام آنالیز حساسیت در 
شبیه‌سازی رواناب رودخانه کرگانرود تعیین شد و همچنین polie‏ حداکثر و حداقل پارامترها به دست آمد. 
در (جدول A‏ پارامترهای محاسبه شده برای مدل‌سازی Gilles‏ با الگوریتم SIMHYD‏ با روش بهینه‌سازی 
جستجوی الگو (pattern Search)‏ نشان داده شده است که بیانگر ٩‏ پارامتر برای تخمین رواناب است. 
ضریب جریان پایه توسط ضرب آن در میزان رواناب سهم جریان اصلی تولید می‌شود. اختلاف Ole‏ بارش و 
تبخیر به‌وسیله ظرفیت ذخیره بارش در برگاب محاسبه شده و سهم نفوذپذیری مقداری از میزان کل بارش 
قابل‌نفوذ در خاک است. لازم به ذکر است polis‏ بهینه‌ای که نزدیک به عدد صفر باشد نشان‌دهندۀ حساسیت 
پایین‌تر پارامترها در bagda‏ موردمطالعه است. بعد از تعیین polis‏ پارامترها. ابتدا داده‌های مربوط به 
سال‌های ۲۰۰۳ تا ۲۰۱۵ برای مرحله واسنجی مدل در نظر گرفته شد و سپس فرآیند کالیبراسیون صورت 
گرفت. برای مقایسه گرافیکی نمودارهای مربوطه برای مرحله واسنجی و اعتبار سنجی ترسیم گردید (شکل ٩‏ 
تا ۱۰). همچنین مقادیر حاصل از نتایج شبیه‌سازی انجام شده توسط مدل SIMHYD‏ در جدول زیر نمایش 


داده شده است (حدول 4 و ۰ A)‏ 


Silas‏ اول 


جغرافیا و مخاطرات محیطی 


YA 


حدول #۸ مقادیر پارامترهای استفاده شده در مدل SIMHYD‏ 
نام پارامتر توضیحات پارامتر حداکثر polis ae‏ بهینه 

۱۸ ۰ ۱ ah Ob > ضریب‎ Base flow coefficient 
*/*0^ ۰ ۵ آستانه نفوذیذیری‎ Impervious Threshold 
۱/۳۷ ۰ tes ضریب نفوذ‎ Infiltration coefficient 
£/AW ۰ E شکل نفوذ‎ Infiltration shape 
۸1 ۰ ۱ ضریب جریان ورودی‎ Interflow coefficient 
۰/۸۵ ۰ ۱ سهم نفوذیذیر‎ Pervious Fraction 
1۰۳۹ ۰ ۵ ظرفیت ذخیره برگاب بارش‎ Rainfall interception store capacity 
UAE ۰ ۱ ضریب تغذیه‎ Recharge coefficient 

Ove \ Ove ظرفیت ذخیره رطوبت خاک‎ Soil moisture store capacity 
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شکا 4- مقادیر ماهانه دیی مشاهداتی 3 شبیه‌سازی شده در مرحله واسنحی مدل SIMHYD‏ 


شکل EE‏ نمودار رگرسیون خطی ميان ee‏ مشاهداتی 3 شبیه‌سازی شده- واسنحی مدل SIMHYD‏ 
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جدول A‏ - مقادیر به دست آمده در فرآیند واسنجی و اعتبارسنجی مدل SIMHYD‏ 


ضریب تبیین(*18) 


ضریب نش اکا 


دوره زمانی 


*/A0NY 


۰/۰ 


واسنجی (۲۰۰7-۲۰۱۵) 


ZEY 


۵ 


اعتبارسنجی (۲۰۱۲-۲۰۱۸) 
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دول ۱۰- مقادیر به دست آمده از ريشه میانگین مربعات خطا در مرحله واسنجی و اعتبارسنحی SIMHYD‏ 


R(XY) | MAPE | RMSE | MSE | MAD n پارامترهای ريشه میانگین مربعات خطا‎ 
۰/۹۲ | ۱ yyy VEY | VAYY | ° مقادیر پارامترها واسنجی‎ 
wav ۱۵/۳۹۲ | WIV | PAY | aiv | Y| اعتبارسنجی‎ 


om‏ برای انجام صحت سنجی مدل شبیه‌سازی شده» از داده‌های مشاهداتی مستفل از مرحله واسنجی 
یعنی مربوط به سال‌های ۲۰۱۲ تا ۲۰۱۸ با همان پارامترهای مؤثر و مقادیر تعیین شده در فرآیند اعتبارسنجی 


در نظر گرفته شد و نتایج به دست آمده به‌صورت نمودار ترسیم شده است (شکل ۱۱ و AVY‏ 
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شکل ۱۱- نمودار دبی مشاهداتی و شبیه‌سازی شده — اعتبارسنجی مدل SIMHYD‏ 


Y = 1/1932 + 77 
R è= 0/9496 
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SIMHYD شکل ۱۲- نمودار رگرسیون خطی میان دبی مشاهداتی و شبیه‌سازی شده- اعتبارسنحی مدل‎ 
SIMHYD نتایج حاصل از مدل‌سازی‎ -۵-۳ 


با توجه به polio‏ به دست آمده از شاحص‌های Goll‏ داده‌های مشاهداتی در واقعیت با داده‌های 
شبیه‌سازی شده و با نگاه اجمالی به نمودارهاء علاوه بر تطابق قابل‌قبول داده‌ها aba‏ پیک و فرود روانات 
بەخوبى شبیه‌سازی شده و ade‏ به دست آمده توسط تابع هدف برای شبیه‌سازی صورت گرفته قابل‌پذیرش 


است. به‌غبارتی‌دیگر با اف ای میزان دبی مشاهداتن عقادیر شبه‌سازی شده نیو افوایش ails‏ و بلعکس با 
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dio as ارف ای ال‎ ole, pod alas, RAS US ceo ا‎ Se کی‎ AS 
در این حوضه. جهت برآورد دبی جریان در تبیین سیلاب نیز از نتایج قابل‎ SIMHYD سنجش کارایی مدل‎ 
مشاهداتی با مقادیر شبیه‌سازی شده حاکی از مناسب بودن‎ polie قبولی برخوردار است. تطبیق تصویری‎ 
مقدار ضریب نش-ساتکلیف(۳) به دست آمده است.‎ 
دبی‎ Ole با توجه به خروجی‌های مراحل واسنجی و اعتبارسنجی و با عنایت به نمودارهاء به‌طورکلی‎ 
مشاهداتی با دبی برآورد شده توسط مدل, همبستگی و تطابق مناسبی وجود دارد و می‌توان این گونه نتیجه‎ 
برای شبیه‌سازی دقیق روند‎ SIMHYD گرفت که با مقایسه نتایج هیدروگراف‌های نمونه باوجود تلاش مدل‎ 
کرگانرود قابل‌قبول‎ ppl همچنان توانایی آن در شبیه‌سازی دبی ماهانه حوضه‎ bw هیدروگراف‌ها در برخی‎ 
است؛ زیرا در شبیه‌سازی ماهانه به دلیل تعدیل در حداکثر و حداقل داده‌ها و نیز از بین رفتن داده‌های نامعتبر‎ 
بالاتری صورت گرفته و نتایج قابل‌قبول‌تری را محاسبه می‌کند. این‎ AS مدل‌سازی به‌درستی و با دقت و‎ 
بسیار‎ pul در ابتدا با داده‌های ورودی روزانه شبیه‌سازی شد و‎ SIMHYD مسئله در حالی اشت که مدل‎ 
ضعیفی را ارزیابی نمود. ازجمله دلایل این موضوع نیز می‌توان به کمبود ایستگاه‌های موردنیاز در حوضة‎ 
موردمطالعه و نیز عدم کیفیت داده‌های روزانة ثبت شده به سبب تازه تأسیس بودن ایستگاه موجود در‎ 
در برآورد رواناب روزانه به جهت نیاز به ورودی‌های‎ SIMHYD کاهش دقت مدل‎ Of محدوده و به دنبال‎ 


کم‌تر نسبت به مدل SWAT‏ ,| اشاره نمود. 
TY‏ ارزیابی نتایج مدل SWAT‏ و مدل SIMHYD‏ در برآورد سیلاب حوضه آبریز کر گانرود 


با توجه به polis‏ به دست آمده توسط مدل‌ها در ماه‌های آخر سال ۲۰۱۱ و همچنین ۲۰۱۶ در مراحل 
واسنجی با افزایش میزان بارندگی. نقاطی به‌عنوان پیک دبی شناسایی شده است که با بررسی حوادث سیل 
حوضه آبریز در این بازه زمانی» مشخص گردید در شهریور و مهرماه سال‌های ۹۰و ٩۳‏ مطابق با پیک‌های 
دبی نمایش داده شده در نمودارهای خروجی. سیل در پایین‌دست حوضه که رودخانه از وسط شهرستان 
تالش می‌گذرد رخ داده و باعث تخریب منازل روستایی واقع در حریم رودخانه. آب‌گرفتگی معابر و CUE‏ 
جانی و مالی بسیاری شده است. در مرحله اعتبارسنجی هر دو مدل نیز به‌ویژه در مدل SWAT‏ که شبیه‌سازی 
دقیق تری را انجام داده است. در سال ۲۰۱۸ نیز با افزایش بارندگی‌های ثبت شده و دبی مشاهداتی» نقاط پیک 
دبی شبیه‌سازی شده است که این امر با سیل ۵ مهرماه سال AV‏ مطابقت می‌کند؛ بنابراین می‌توان نتیجه گرفت 
polis‏ به دست آمده توسط مدل SWAT‏ و مدل SIMHYD‏ در برآورد حجم رواناب ایجاد شده به‌حوبی 


عمل کرده و نقاط سیل را به‌خوبی شبیه‌سازی نموده‌اند (شکل ۱۳). 
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شکل ۱۲- مقایسه نتایج شبیه‌سازی ماهانه مدل‌های SWAT‏ و SIMHYD‏ با دبی مشاهداتی 


é‏ جمع‌بندی 
با توجه به تنوع مدل‌های بارش -رواناب انتخاب یک مدل مناسب برای حوضه‌های آبریز جهت افزایش 
بهره‌وری برنامه‌ریزی‌ها و مدیریت منابع آبی» امری مهم و ضروری است. به همین جهت در تحقیق حاضر نیز 
برای شبیه‌سازی سیلاب از دو مدل SWAT‏ و SIMHYD‏ بهره گرفته شده که با توجه به نتایج به دست آمده 
تطبیق مدل‌ها با حوضة موردمطالعه از خروجی‌های قابل‌قبول و رضایت‌بخشی به شکل ماهانه برخوردار 
است. به‌طورکلی هماهنگی و تطبیق نتایج مدل‌ها با شرایط واقعی حوضه آبریز کرگانرود. با وجود نواقص 
داده‌های مشاهداتی» مقادیر را تا حدودی به‌درستی اندازه‌گیری کرده و می‌توان برای اقدامات مدیریتی جهت 
تشخیص برآورد و مهار سیلاب با توجه به رواناب تولید شده به نتایج این مدل‌ها اتکا نمود. QS‏ و 
همکاران(ع۱۳۹). محضری و همکاران(۱۳۹۵). ابراهیمی و همکاران(۱۳۹۷) سفیان و ماگاجی" (۲۰۱۹) و 
فیری و همکاران" (۲۰۲۱) نیز در حوضه‌هایی مشابه در شمال کشور با استفاده از مدل SWAT‏ رواناب را 
شبیه‌سازی نموده و نتایجی مشابه پژوهش انجام گرفته به دست آورده‌اند. همچنین چیو و همکاران " (۲۰۱۸)؛ 
موبیالیوو و همکاران" (۲۰۲۱) و gd‏ و همکاران" (۲۰۲۲) نیز با استفاده از مدل(۹1۷1۳17۳ در بررسی سیل به 
نتایج خوبی دست یافتند. در این میان لازم به ذکر است که برای این پژوهش مدل SWAT‏ به‌عنوان ابزاری 
قدرتمند. عملکرد نسبتاً بهتری در مقایسه با مدل توده‌ای SIMHYD‏ در ارزیابی تولید رواناب ایجاد شده در 
این حوضه داشته است. ازجمله دلایل این برتری را می‌توان ماهیت نیمه‌توزیعی بودن آن و نیاز مدل SWAT‏ 


به ورودی‌های بیش‌تری نسبت به مدل SIMHYD‏ ذکر کرد که منجر به شبیه‌سازی رواناب روزانۀ این حوضه 


1 Sufiyan and Magaji 
2 Phiri et al 

3 Chiew et al 

4 Mubialiwo et al 

5 Liu et al 
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ترخات Urey Gaol utd SIMY Da‏ عامل مهس os‏ ایفاه یلیهان جو فة x‏ € 33318 
بشمار می‌رود. لذا به سبب این ماهیت L‏ عدم استفاده از مقادیر ثابت در ورودی‌هاء دقت مدل در شبیه‌سازی 
حوضه‌آبریز مطابق با واقعیت افزایش می‌یابد. از طرفی دیگر بر اساس نتایج به دست آمده توانایی یک مدل 
در ارزیابی یک حوضه آبریز تنها در میزان پیچیدگی Ol‏ خلاصه نمی‌شود بلکه کیفیت و اعتبار داده‌های 
ورودی مدل در درجه اول حائز اهمیت است و خروجی‌های به دست آمده در هر دو مدل بیانگر این موضوع 
می‌باشند. چراکه حوضه‌های آبریز کشورهای دیگر به‌ویژه امریکا و استرالیا در حوضه‌های مشابه توسط همین 


مدل‌ها برآورد شده است و نتایج بسیار خوبی را مطابق با واقعیت شبیه‌سازی نموده‌اند. 
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